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Abstract 



Revivals of optical coherence of molecular photoassociation driven by two 
ultrashort laser pulses are addressed in the Condon approach. Based on textbook 
examples and numerical simulation of KrF excimer molecules, a prediction is 
made about an existence of photon echo on free-bound transitions. Delayed rise 
and fall of nonlinear polarization in the half-collisions are to be resulted from the 
resonant quantum states interference whether it be in gas, liquid or solid phases. 
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1. Introduction 

At present, the photoassociation spectroscopy (PAS) of atoms colliding along 
the ground state potential and interacting with photons to produce molecular 
dimers (through the inverse predissociation) has attracted a great deal of at- 
tentions This technique, as applied in conventional continuous wave^'^^ 
(CW) and fentosecond^ regimes, permits a determination of the rovibrational 
leveles, hyperfine structure, scattering lengths and/or dynamics of bond forma- 
tion that can be derived with an unprecedented accuracy. To assess the primary 
processes of harpooning reactions^ , photoassociative and the Penning ioniza- 
tions ^, relaxation of the van der Waals complexes^ in liquids and quantum 
gases in optical traps, the PAS is ideally suited. The review of earher devel- 
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opments, which have been collected under the rubric of chemical reactions in 
transition spectroscopy, can be found in ^. The most recent innovation to the 
CW PAS is investigation of Bose's state of the cold atoms^ in which a loss of 
coherence can be controlled trough the laser assisted binary collisions. 

One would naively claim that the femtosecond PAS is not quite sensitive 
as the resonant CW schemes, since only a low number of atoms distributed 
randomly in close proximity during ultrashort pulses may contribute to the co- 
herent bonding. However, the shorter impulse, the wider its spectrum, the more 
the free-bound transitions involved. An example is found in photoinduced har- 
pooning reactions ^ in which according to the M. Polanyi the attacking open 
shell atom tosses out its valence electron, then hooks the closed shell atom and 
hauls it in with the Coulomb force even at a range of hundreds Angstroms^. 
This reaction endowed with separation of degrees of freedom into nuclear and 
electronic coordinates controlled by the femtosecond optical pulses may have 
a high quantum yield. 

It is straightforward to extend the femtosecond PAS to address the four wave 
mixing geometry or its two or three wave realizations. This makes it possible 
to investigate the nonUnear optical effects those are no longer observable in the 
weak fields. Such as a photon echo photoassociation which is counterpart of 
the photon echo of photodissociation. A theoretical description of the latter is 
only being reported by us ^'^^ in context of the wave packet engineering. It 
is evident now, that a well matched spreading and squeezing of the vibrational 
wave packets are called for the desired optical properties of molecular samples. 
So that the quantum control in the weak field regime requires a complicated 
pulse shaping While "spontaneous" coherence, i.e., a photon echo, is ob- 
tained trough a more strong, short and transform limited (non chirped) fields 
that enables one to automate the wave packet preparation on the free-bound 
transitions. 

The transients having the entanglement of discrete levels and continuum 
usually result in a decay The quantum transition ampUtudes may decrease 
depending on the contribution of the molecular repulsion or quantum disper- 
sion. In the frequency domain, the resonant coupling is accompanied by the 
inhomogeneous broadening (dephasing) of the Franck-Condon transitions. If 
the latter dominates over an irreversible relaxation in multi-channel (chemical) 
collisions, spontaneous emission, the photon echo can be harvested as for the 
bound-bound transitions 

To make this point a bit clear, we follow the next plan. In Sec.2 the Condon 
model of femtosecond PAS is analyzed. We relate the delayed transients to the 
wave packet interference residing in molecular continuum after the second im- 
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pulse. Analytical expressions are derived for the x^^^ polarization induced by 
2 (5-like impulses. A singularity of x*^^^ is identified with an infinite region of 
the Franck-Condon transitions simultaneously involved. For the realistic field 
shapes this trait is eliminated. The PAS is particular appealing for the investi- 
gating excitations of the rare gas halides, for which the electronic ground states 
are dissociating. This system is of use as an active media of the high-power 
lasers^ and photon echo of photoassociation might deliver a valuable informa- 
tion on timescales of excited molecular states in the lasing. Specifically, the 
relevant model of molecular dynamics of the krypton fluoride KrF(B) excimer 
is simulated in Sec. 3, where we consider a photoinduced harpooning reaction 
of interacting Kr cations and F anions created at thermal conditions. To avoid 
the Condon reflection principle, which overestimates contribution of turning 
points and disregards oscillatory tails of the wave functions, the B-X transition 
(at 248 nm) is treated in the globally uniform approach. The Wigner represen- 
tation gives an insight into why the photon echoes of photoassociation exist 
from the kinematical point of view. Section 4 concludes. 

2. Condon Model 

Strictly speaking, 'ab-initio' approach to PAS necessitates a many-body 
field theory of interacting matter-light species. However, any more rigorous 
treatment of the photoassisted reactions beyond the Born-Oppengheimer strat- 
egy faces sophisticated mathematical problems of exclusion of extra degrees 
of freedoms and transitions between multidimensional surfaces. It is not sur- 
prising, then, that a generic model with simplified potential curves, in which 
the bound-free transitions occur in one-body treatment, has been introduced at 
the drawn of the age of quantum mechanics Herein, a separation of degrees 
of freedom gives the paradigm of application of the Franck-Condon principle 
asserting that for a time of electron transfer, the coordinates and momenta of 
nuclei do not vary because of the large difference between the electronic and 
nuclear masses. Molecular dynamics is governed by the Hamiltonian opera- 
tors: 

Hb = Tkin + Vb', Hf = Tkin + Vf, 

where T^in is the operator of kinetic energy of relative motion in the motionless 
center mass of a molecule and the Doppler's shift of frequency is disregarded 
as marginal. Initially, atoms collide along the ground electronic term with 
constant interatomic interaction Vf = 0. The excited electronic curve can be 
presented by the potential Vf, = ujo+Q'^x'^ /2 with understanding that the actual 
curves are deflected of the parabolic well at longer and shorter distances from 
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its equilibrium position. The model potential has merit of being exactly solved 
best serving for our methodological purpose. We have chosen the test example 
of oscillatory motion with frequency 17 as a vehicle with which to demonstrate 
general aspects of photoassociation. Without restricting the generality, an ef- 
fective mass of molecule equals unit. The electronic states are coupled by the 
electric dipole interaction Vbj = fiE{t), where the laser impulses have the 
slow envelops E{t) and fast optical frequency ujq. The envelope differs from 
zero only for the moments of photon impacts. 

With the rotating wave approximation, the Schredinger equation for the 
bound states IV'fe) and scattering states can be written in the notation of 
Dirac as 

in\'iPb) = Hb\i>b) + Vbf\i>f) 

in\i^f) = i?/|V/) + v>6|V6)- (1) 

In order to have general formalism capable of dealing with both bound and 
scattering states, we employ the momentum representation. A crucial assump- 
tion is made that our system has a very large box volume V. The energy spec- 
trum of atoms forming the plane waves is bound to be quasi-continuous in the 
box. It is required of the initial quantum conditions (^1^^/(0)) = V'j(p) that 
the box wave functions possess the property of random phases. This require- 
ment corresponds to the "molecular chaos" with the canonical Gibs thermal 
distribution in the atomic ensemble: 



Pi J 



where (3 is the inverse temperature, Z is the statistical sum. The box normaliza- 
tion (a molecule per the box volume) is implied. The quasi-continuous states 
correlate initially only with the coincident momenta giving the Gibbs density, 
while the bound molecular states are empty |V'6(0)) = before the laser exci- 
tation. 

The resonant states between the laser impulses are governed by the molec- 
ular evolution operators 

\Mt)) = e-*^//'^|V'/(0)), (2) 

where the momentum representation can be used to arrive at the matrix ele- 
ments 

. tii^ gi((p2+p2)cos(m)-2ppi)/(2?insm(m)) 
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The laser field couples the resonant electronic states IV'fe) and \il)f) and to pre- 
pare their mixing in the Umit of a weak, short impulse according to eq. (1) 

Hf{r)) = H))-i9mr)). (4) 

where zero time t = corresponds to the beginning of an ultrashort pulse of 
the duration r, so that 9 = Tjj, ■ E jU <^ it hihe pulse area. 

Let us consider an idealization of the 5- like impulses. The Franck-Condon 
principle ensures the intact positions and momenta of the wave packets on the 
resonant electronic terms during the laser field action. The optical polariza- 
tion of the system is found by calculating the dipole moment pL in the Condon 
approximation in which a total span of Franck-Condon region of free-bound 
transition includes the (infinite) box V, whilst the relative electronic transition 
moment, i.e., the dipole operator jl is independent of a position of atoms in the 
box: 

{m)ave = n{{lpf{t)\Mt)))ave- (5) 

The free decay of the optical polarization after the first laser impulse follows 
the time dependent Fermi golden rule: 

(A^^^)ar-e = -iiieTr[e-'^^»/h^i{t+if3)Hf/hyz = -i^0K(^){t), (6) 

where the indentity relation / = / dp\p) {p\ and the matrix elements of Eq. (3) 
are appUed for the trace calculation 

= Z-^[cOs{Wdecay{t))-l]-^/^. (7) 

The kernel K^^^ is determined by a "hyperbolic" rotation on complex angle 
W{t) given by: 

cos{Wdecay{t)) = cos{nt) + 0.5n{i/3 + t)sin{nt), (8) 

The hyperbolic rotation is inspired by the dynamical symmetry group SU{1, 1) 
responsible for the phase and ampUtude modulation of average dipole mo- 
ment. In frequency domain, the hyperbolic rotation gives rise in inhomoge- 
neous broadening of spectral lines of free-bound transitions 

The free-bound transition amplitude K^^^ (t) has a time-integrable singular- 
ity of the order 1 / s/t. Notice that this kernel is resulted from an unbounded 
spectrum of the simultaneously driven resonant transitions induced by the 6- 
pulse. This artifact of the polarization singularity can be corrected whether 
by a dipole moment cutoff that leads us beyond the Condon approximation 



6 



and/or frequency limited spectrum of the excitation field. However, being the 
/i cutoff length dependent, it becomes tedious for the analytical analysis. The 
thermalization and internal molecular interference add to the complexity of the 
problem. This is the case of femtosecond PAS. The nonlinear optical polar- 
ization depends also on the multi-dimensional convolutions of the actual field 
envelopes with the singular kernels induced by the 5 fields. Their nonlinear 
combinations provide another vectors of directional coherence along which 
the echo signals are to be observed. We shall set these problems aside for 
sake of simplicity to treat molecular photoassociation in the (5-pulse limit and 
focusing on collinear geometry of incident electromagnetic wave vectors. Sug- 
gesting that these two 5-pulses operate having the inter-pulse delay T, we find 
the mixing bound and free states as 

\MT)) = -ie^e-'^^'f/^l^j), (9) 

where the initial state |V'°) belongs to the Gibbs distribution of collision pairs. 

By 9i and 02 we denote the areas of the first and the second laser impulses 
respectively. These pulses govern a quantum superposition of the states 

= -i02\Mt)), (10) 

in which the normalized state vectors \'ipb{t)), \'^f{t)) entering Eq.(lO) with 
the coefficients proportional to 02 and 9i92 yield a dominant contribution to 
the delayed transition amplitudes. We call these vectors as the "propagating" 
states The free-bound transients are presented by a series of the propaga- 
tion operators each of them being responsible for an adiabatic step along the 
molecular curves. Therefore, a first nonlinear correction to the dipole moment 
averaged over the Gibbs distribution can be cast as 

{fi^^^)ave = if^9i9jK^^\t), 

/o\ ■(T-t)Ht, .(t-T)Hf ,THf, (g-iT)Hf 

= rr[e*— ^e*— ^e-*-^e ^]/Z. (11) 

Using the momentum representation of the propagation operators and perform- 
ing this 2-d Gauss integration over momenta in the matrix elements Eqs.(3,13), 
we obtain a characteristic function 
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Figure 1 The free radiation 
decay and echo transients in 
the Condon model. The 
higher temperature, the faster 
decay and the sharper the pho- 
ton echo peaks. This be- 
haviour is illustrated by two 
curves at inverse temperature 
p = 1/n and P = l/5fi. 
Inset displays the //^^^ dipole 
moment arising on free-bound 
transition following the sec- 
ond ultrashot impulse which 
fires with delay T « 2n/il 
matching molecular period. 
Use of normalized impulse ar- 
eas 01,2 = 1 and frequency 
f2 = 1 is made here. 

in which a hyperbohc rotation owing to the 2-pulse interaction is 

cos{Wecho{t)) = cos{nt} + 0.5n{{i + il3)sin{nt} 

+ n{i+T){T -ip)sin{n(t + T))sin{nT)), (13) 

being conveniently expressed trough the shifted time t = t — 2T. 

For the weak field regime the resonant coupling of continuum and bound 
wave functions in the Franck-Condon region is determined by the time de- 
pendent Fermi golden rule as shown in Fig. 1. A monotonic decay of the 
Franck-Condon factors will no longer be valid at a low temperature, e.g. at 
inverse temperature /3 = 1 /fi corresponding to one vibrational quantum. The 
reason has to do with the quantum dispersion which alone cannot destroy co- 
herence. Molecular repulsion for inclined potential curves or predissociation 
would afford exponential decay But even in this instance, a net construc- 
tive or destructive interference of two wave packets prepared by the pump and 
delayed impulses occurring at or near multiples of the vibrational period might 
result in the revivals of coherence. The transition amplitude spikes reflect the 
return of the vibrational wave packet on the attractive potential curve. 

To consider the echo effect in the weak field regime we derived the 
polarization from Eqs.(l 1-13). In the Condon model, the transition amplitudes 
between the propagating states are obtained explicitly as a function of time 
and an average polarization is plotted in inset in Fig.l. The echo signal itself 
consists of a number of short peaks most prominent when a delay T is close to 
multiple of the vibrational period (T f» = 2it/Q). The Eq.(13) determines 
the positions of the "spontaneous" revivals near t = nT, (n + 1/2)T n = 
1, 2.... In the occasion \sin{ilT)\ <C 1 and (Orsm(Or))^ > 1, we easily see 
that the half width at a half maximum (HWHM) att = 2T is estimated from 



Photon echoes of molecular photoassociation 




5 r >. h S II) 12 14 1h tg 



rinic |ix'L uiiitsj 




1 2 3 4 5 6 



Time [rel mniti] 



8 



Eq.(13) as 

Techo ~ (14) 

This simple expression can be understood as a time taken by a classical par- 
ticle to move the thermal de-Broglie wavelength rp oc with velocity v = 
Q.{R)ave depending directly on the correlation radius {R)ave of the propagat- 
ing state in the continuum of states, 

{R{2T))ave = {{i^f{2T)\R\^f{2T)))ave OC Tvo, 

where vq ccl/ vT? is the thermal velocity of atoms. So that, we obtain the cor- 
relation time Tcor = ro/TvQ = P/0,T, which matches Eq.(14) for HWHM and 
has to play a key role in determining the nonlinear polarization peaks. What 
this means is the validity of kinematical argumentation based on the quasiclas- 
sical transport theory allowing us to justify the echo of photoassociation just 
as the photon echo of molecular photodissociation^'^ . It is appropriate to re- 
call that in the latter, the ground molecular wave function forms a wave packet 
of the receding atoms following first impulse. Atomic distribution acquires a 
stretched elliptic shape in the phase space formed by the coordinates and ve- 
locities. The ellipse is strongly squeezed in the transversal direction by virtue 
of "free" propagation. However, the total phase space volume (according to 
the Liouville theorem) persists. Being driven back by the second impulse, the 
wave packet replica settles in the binding term as a non stationary state. 

For the harmonic model, the ellipse rotates as a hole at the vibrational fre- 
quency. In addition, the replica of the initial state, which was preserved after 
the first impulse, repeats the evolution of the wave packet in the molecular con- 
tinuum. The transition amplitude \K^^^ p strongly depends on the interference 
between the quantum-mechanical paths of the wave packets driven in tandem. 
The last manifests in squeezing of the photon echo peaks inversely quadratic 
of the delay T: 

Techo = Tdecay/mTf + 1), (15) 

where Tdecay is the time of free radiation decay depending on molecular curves 
and temperature. It equals a molecular radius a divided by typical velocity 
of receding atoms in the phase space. The longer the inter-pulse delay T, the 
larger the squeezing of molecular states, the larger the correlation radius of 
receding atoms, and the faster decay of the Franck-Condon overlapping. This 
implies, that the wave packet moving on the smaller ellipse diameter a varying 
as 1/ QT has to move it with velocity v which grows directly QT (proportional 
to the correlation radius). The overlap time is given by the ratio a/v to yield 
the mentioned trend Techo oc (ilT)~^. 

The same reason holds true for the photon echo of photoassociation. But its 
duration features its inverse power of T. This is because only a propagating 
state in the continuum is being squeezed after the second impulse as a function 
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of inter-pulse delay. The difference between the echo scaling laws for the pho- 
todissociation and photoassociation is founded on the fact, that, in the former, 
both wave packets evolved on the binding term and in the continuum X-states 
are squeezed. Also the ji^^^- polarization has a maximum peak near t = 2T. 
Our analysis indicates that the HWHM increasingly decreases with the temper- 
ature. Besides, at the delay time exactly matching the oscillatory period 27r/Q, 
the amplitude of spontaneous revivals may appear to diverge at points t = 2T. 
Fortunately, this divergence rules out in the /j,^^^ polarization controlled by the 
realistic laser field which shape differs from the (5-like envelopes. 

The theory of this phenomenon is best clarified in the phase space of the 
system. The free-bound transitions are characterized by the distribution of the 
coUisional pairs over their electronic excitations, intemuclear separations and 
velocities. Therefore an uniform representation for the scattering states, those 
are quasiclassical in their nature, may be helpful. Just as for the photodis- 
sociation picture, the molecular dynamics can be visualized for the Wigner 
distribution functions. For the thermal atoms, the Wigner map has the shape 
of infinite strip along the coordinate axis. The strip width in the transversal 
direction is equal to the thermal de-Broglie wavelength. 

The first impulse conveys the coUisional pair to the binding electronic B- 
term. The replica of the thermal Wigner distribution sets here. After that, the 
atoms oscillate between the turning points. This corresponds to the strip dis- 
tribution in phase space which rotates at the harmonic frequency. The second 
impulse does the same with the remaining portion of colliding atoms. In ad- 
dition, it simultaneously promotes the rotated strip to the continuum, where 
the latter starts to spread along the coordinate axis. The spreading distribution 
squeezes in transversal direction. It is apparent, that their overlap sharply rises 
and falls, when the adjacent Wigner maps cross each other. If the maps turn out 
parallel, a total span of the induced free-bound transitions includes an infinite 
volume of our simplistic model enabling the Franck-Condon factor to diverge. 

3. Wigner Functions of Photoassociation of Kr-F 
CoUisional Pairs 

Consider the product of the density matrices of the ground and excited elec- 
tronic states. The time-delayed nonvanishing overlap of the states in phase 
space would signature an existence of photon echo of photoassociation. For 
ultrashort pulses, the propagating states contribute to the delayed transients. 
Then, the density matrices pb = {\tpb) {ipb\) ave and pf = {\ipf){ipf\)ave are 
related to the Wigner distribution function Wa{p, q) via the Fourier transform 
in the coordinate space 



Wa{p,q) = {2Tr)-^ j dre'P'{{q-r/2\iPa) O'alq + r/2))ave, a = b,f.(\6) 
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Figure 2 Interaction be- 
tween Rr cation and F anion 

is described by the truncated 
Ritner potential^®, Vb{r) = 
a + fee — cir^^ — C4r^'^, 
where r is the Kr-F dis- 
tance, a = 85510 [cm~^] 
is the separated ion limit, 
7 = 2.94 [A-^] and 
b = 1.01 X lO'' [cm-^] 
are the constants of electron 
shells exchange repulsion, 
ci = 2.19X IQ-'"' [cm^M] and 
C4 = 1.7 X 10^ [cm-'^A*] are 
the parameters of Coulom- 
bic and ion-quadmpole 
interactions respectively for 
the SS+,2 state. For the 
ground -'fS+^j state, the 
interaction potential is as- 
sumed to be purely repulsive: 
Vx{r) = cisr"^^, where 
ci3 =3.2 X 10^[cm-^A^^]. 



The overlap of the Wigner distribution functions in the phase-space is ex- 
pressed as 



Tr[pbPf] 



dpdq Wb{p,q)Wf{p,q). 



(17) 



In our simple model, the quasiclassical functions Wa {p, q) follow the Liouville 
transport equations 



+ P 



dVa dWa 



dt 



dq dq dp 



= 0, 



(18) 



with an accuracy of independent of the model potentials under study. This 
approach guarantees an uniform treatment of molecular dynamics even though 
a major contribution to the overlap integral is accumulated near the turning 
points, where the standard quasiclassical theory fails. We take an advantage of 
the Wigner representation to calculate the matrix density overlap of Kr(4p^S'o)- 
F(2p^ ^P) photoassisted collisions along X- and B -molecular terms shown in 
Fig.2. Applying the classical dynamic routines to colliding Kr and F atoms, 
we directly obtain the Wigner density functions of photoassociation. The 5- 
like envelopes of laser fields allow us to deal just with the propagating states 
density. In Fig. 3 we display the contour maps calculated from the transport 
equation Eq. (17). The thermal atoms occupy the half-strip terminating in a 
short distance due to the strong exchange repulsion on the covalent X-term. 
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Figure 3. The Wigner contour map of the density of KrF excimer molecules immediately 
following an ultrashort pulse at fs, 50 fs and 90 fs. Upper panels show evolution of the Wigner 
functions due to a harpooning reaction on the ionic B-curve. Lower panels show the thermal 
distribution of the colliding atoms along the covalent X-curve. The contour lines are taken at 
half maximum height of the Wigner's functions 
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The energy gap between the ground covalent state and bottom of the ionic B- 
term is in range of 248 nm wavelength resonant to the laser pulse frequency. 
The Wigner maps in the upper row in Fig. 3 present the evolution of the thermal 
distribution of the X-states distorted by the first 5-pulse that resonantly drives 
them to ionic B-states. Then, the F anions and Kr cations are attracted together 
due to the joint action of the Coulomb force and induced ion-quadrupole in- 
teractions. In the harpooning reaction, the Wigner map of the thermal atoms 
curves forming a hook and then a spiral structure in the phase space. The spi- 
ral step is conditioned by the potential anharmonicity. For the pure harmonic 
model, a rotation of the initial strip as a hole at the oscillatory frequency takes 
place. The Wigner maps of the propagating X,B-states created by the pulse 
pair are shown in Fig. 4. We take the inter-pulse delay equal to about a period 
( 90 fs) of molecular vibration in the bottom of B-state. The second impulse 
transports the Wigner's B-state distribution partially back to continuum, as it 
is, conserving its spiral shape. The spiral branches are spreading on the co- 
valent X-curve and simultaneously squeezing in the momentum direction. In 
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Figure 4. The Wigner distribution functions of KrF excimer created by the impulse pair. The 
interpulse delay is equal to 90 fs i.e. about one molecular period. The Wigner maps are drawn 
at 150 fs, 200 fs and 280 fs. Upper panels present the harpooning reaction in the ioiuc B-states. 
Lower panels show the spreading photofragments in covalent X-state. 
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addition, its optical field might induce a haipooning reaction by producing a 
collisional pair of cation Kr and anion F. Tiie Wigner distribution on their B- 
binding term repeats the evolution following the first impulse. In Fig. 5 we 
display the overlap of the molecular densities in the phase space as a func- 
tion of time. The overlap of the Wigner's distributions of the X and B -states 
quickly tends to zero. The decay time after the first impulse can be evaluated 
as a width at the half maximum height of the integral Eq.(17). This time is 
about 20 fs long that conforms the estimate based on the classical mechanics. 
Doing so, consider a particle moving with a thermal velocity in the covalent 
X-state. Its wave packet is transported to the binding B-term via a vertical 
Franck-Condon's transitions. Then, it needs the time: 

/•Rout 

Tdecay ~ / dr/vb{r) (19) 

•' Rin 

to destroy an alignment of the Wigner's distributions while the wave packet 
travels between the turning points in the quantum well. The turning points 
Rin 2.15^ and Rout ~ 2. 65 A are determined by the thermal energy 
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Figure 5. The free radiation decay and photon echo of molecular photoassociation of KrF. 
The HWHM of ijf-^^ polarization is about 20 fs. The relevant /x*-^' polarization created by two 
5-like laser fields having nomialized pulse areas ^1,2 = 1 and inter-pulse delay 90 fs is shown 
in inset. 




Time [fs] 



Ether = 100 [cm~^] of colUsion pairs of Kr and F atoms. The wave packet 
loses its momentum correlation with nascent replica of X state for time its 
journey between the turning points. Our first rough estimate of about 40 
fs is obtained by substituting to Eq.(19) the velocity: vi){R) = {2{Ether — 
Vb(i?))/M)^/^, and by taking the numerical values of the reduced mass M 
15.5 of the KrF molecule. Inset also shows the rise and fall of the overlaps 
of the Wigner's distributions after the second ultrashort pulse. This occurs at 
twice time-delay, when the hook map crosses the spiral branches of the X-state 
distribution. Our numerical calculation indicates, that the photon echo signal is 
rather robust with respect to a small deviations of the inter-pulse delay from the 
gating conditions T = T^if,. This robustness can be attributed to anharmonicity 
of the KrF(B) potential. 

4. Conclussions 

We have demonstrated how the ultrashort pulses can reverse the natural 
tendency of the nonlinear optical transients to decay in photoassisted half- 
collisions. The photon echo of molecular photoassociation presents a paradigm 
of the coherent self-organization in optics. This occurs in the continuum states 
initially residing in the Gibbs equihbrium ensemble. The photoassosiation is 
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followed by the delayed photodissociation giving rise the correlation radius of 
propagating states to grow proportionally as the inter-pulse time-delay. The 
latter results in a long lived transient complex. Then, the echo signals man- 
ifests on the free-bound transitions following the free radiation decay that is 
dictated by the resonant coupling probability in accordance with the time de- 
pendent Fermi golden rule. This delayed coherence is accounted for by the 
inverse predissociation amounting to phenomena conceptually identical with 
the interference stabilization of atoms (molecules) in the strong ionizing (dis- 
sociating) laser fields. 

The femtosecond coherence of the laser assisted collisions is formed on 
the molecular nano-scales and can be observed in optically thin samples. This 
effect might provide a way to make a compact light compressor without limita- 
tions imposed by the self-modulation effects in the optical wave-guides. From 
the spectroscopic point of view, the photon echo of photoassociation provides 
for unique data on the fly during the binary half-collisions. The echo spikes 
splitting can exhibit an inhomogeneous broadening inherent to the distribu- 
tion of guest atoms in host solid with the sub-wavelength resolution up to 
Angstrom's scale. This opens up a possibility to trace the migration of open- 
shell atoms in rare gas matrices and cage exit. Kindred phenomena might be 
studied by molecular dynamics simulations to provide a new avenue for the 
echo spectroscopy in solids and liquids. 

In brief, therein our major aim was not a quantitative agreement with an 
experiment, but to broadly predict the general aspects of the photon echo of 
photoassociation and to delineate those its features which depend on inter- 
pulse delay of laser pumping, molecular potentials and temperature of sample. 
The echo scheme has a merit to automate the wave packet preparation mak- 
ing it a quite favorable for a practical design, where the pulse shaping is not 
necessary. It is important to stress a geometrical aspect of the effect, in which 
the directional vector ki — 2k2 of the echo coherence is distinguished from the 
incident wave vectors ki and k2 of the pump and delayed fields. This offers a 
convenient experimental condition for observing spontaneous photons on dark 
background with angular and temporal lag from the laser firing. The details of 
the signal may be varied from one molecular system to another, but the broad 
conclusion is quite robust. 
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